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Abstract. Deriving the motion of a compact mass or charge can be compUcated 
by the presence of large self-fields. Simplifications are known to arise when these 
fields are split into two parts in the so-called Detweiler- Whiting decomposition. 
One component satisfies vacuum field equations, while the other does not. The 
force and torque exerted by the (often ignored) inhomogeneous "S-type" portion is 
analyzed here for extended scalar charges in curved spacetimes. If the geometry 
is sufficiently smooth, it is found to introduce effective shifts in all multipole 
moments of the body's stress-energy tensor. This greatly expands the validity 
of statements that the homogeneous R field determines the self-force and self- 
torque up to renormalization effects. The forces and torques exerted by the 
S field directly measure the degree to which a spacetime fails to admit Killing 
vectors inside the body. A number of mathematical results related to the use of 
generalized Killing fields are therefore derived, and may be of wider interest. As 
an example of their application, the effective shift in the quadrupole moment of 
a charge's stress-energy tensor is explicitly computed to lowest nontrivial order. 
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1. Introduction 

Essentially all observational consequences of general relativity rely in some way on 
understanding how objects move in response to the curvature of spacetime. This 
problem has been extensively studied in many contexts. One area in which some 
uncertainties still remain relates to the "bulk" interaction of matter with its own 
(gravitational or other) long-range fields. Self-force effects such as these were first 
studied in flat spacetime electromagnetism. At the time, the main motivation was to 
explain the then newly-discovered electron. Today, interest in the problem has been 
revived in order to better understand the gravitational self-force involved in the decay 
of astrophysical compact binaries. 

The first work on the electromagnetic self-interaction problem involved direct 
calculations of the force acting on a charged rigid sphere in flat spacetime [H El l3| ■ At 
least in the nonrelativistic limit, the acceleration a of the sphere's center was found 
to be closely approximated by 

2 

md — Fcxt + -z'f'o, — Sma. (f ) 
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q represents the body's total charge and m its mass (as calculated by integrating 
over the interior in the usual way). -Fcxt is the ordinary force exerted by the external 
electromagnetic field. The last two terms on the right-hand side of ([T]) are due to 
the self-field. The first represents the expected reaction to emitted radiation. The 
other effectively shifts the charge's inertial mass. The formula for 6m is relatively 
complicated. It has no meaningful point particle limit, and is strongly dependent on 
the details of the charge's internal structure. 

For these reasons and others, it has been common to interpret the mass 
renormalization term as "practically unobservable" in most contexts. The inertial 
mass of a particle is usually inferred by observing its reaction to applied forces. It 
is clear from the above equation that such a mass is given by nicft = m -\- 5m. The 
relation between m^s and the body's internal structure is often irrelevant. It is just 
some parameter that can be fixed experimentally (at least if there is reason to expect 
rhoff = 0). 

This point of view has led to methods that directly extract only the more 
interesting "radiative" portion of the self-force. One proposal - originally due to 
Dirac [4] - is to compute the force exerted by one-half of the difference between 
the charge's retarded and advanced fields. The result satisfies the vacuum Maxwell 
equations, and almost always varies slowly inside the body. It is also insensitive to 
the details of the charge's internal structure. A simple calculation of the force that 
it exerts recovers the classical Abraham-Lorentz-Dirac expression with far less work 
than most other methods. Similar ideas have recently been generalized by Detweiler 
and Whiting in order to recover standard results for scalar, electromagnetic, and 
gravitational self-forces in curved backgrounds [5]. In any of its incarnations, the 
effects of self-interaction on a body's motion can be interpreted as arising only from 
a particular vacuum component of its self- field (up to "renormalization effects" ) . 

We will refer to this idea as the Detweiler- Whiting axiom. It was originally applied 
to point particles as a physical explanation and a shortcut to the conclusions of various 
perturbative calculations. It was not, however, derived in any rigorous sense. Its 
accuracy and range of validity were therefore unknown. This situation has recently 
changed. A direct derivation of a suitably interpreted Detweiler- Whiting axiom is 
now available in very general contexts. Natural linear and angular momenta were 
non-perturbatively defined for nearly arbitrary extended bodies in curved spacetimes. 
These were shown to have evolution equations that depended only weakly on the 
inhomogeneous "S-type" self-field [51 [71 [5] . The main effect of this field was absorbed 
into the definitions of the momenta; a direct generalization of the mass renormalization 
effect illustrated by IJ) . It is the only consequence of the S field in Minkowski and de 
Sitter backgrounds. More generally, the existence of any single Killing vector implies 
that the evolution equation for a particular component of the momentum is completely 
independent of the S field. If there are no background symmetries, a similar result can 
still be derived to lowest order in a standard approximation scheme. The S component 
of the self-field has a small effect in this case, however. It does more than shift the 
effective momenta. 

This paper discusses the nature of these extra forces. Mathematically, they arise 
as a measurement of the degree to which a particular Green function fails to be 
invariant under the action of certain generalized Killing fields. The result of this 
broken symmetry is shown to have a very simple interpretation. The quadrupole and 
higher multipole moments of a body's stress-energy tensor are effectively shifted by 
small amounts. In combination with the aforementioned momentum renormalization. 
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this demonstrates that the only direct effect of the S field is to make an object move 
as though it had a different stress-energy tensor. The effective stress-energy tensor 
couples to the gravitational tidal fields, just as it would even for a test body in 
Newtonian theory. The self-field therefore infiuences a body's "mass distribution" 
in both inertial and (at least passively) gravitational senses. 

The framework used to derive these results has been developed for objects 
individually coupled to scalar, electromagnetic, or (linearized) gravitational fields 
[H El [8]. A set of generalized Killing vectors must be defined first 9 . This 
requires a preferred worldline and a set of spacelike hypersurfaces foliating the 
object's worldtube. Both of these constructions are usually fixed using center of mass 
conditions. Linear and angular momenta for general compact matter distributions are 
then given as conjugate to each generalized Killing field. Their evolution equations 
are determined using stress-energy conservation. If the homogeneous R-type self-field 
is not too large, the laws of motion in this formalism coincide with those derived by 
Dixon [inilllllll] in the test body hmit. 

These ideas are not widely known, so a review of the theory as it applies to 
scalar fields is presented in Sect. [21 Those results are then used in Sect. |3lto argue 
that a charge's S field only acts to (finitely) renormalize the multipole moments of 
its stress-energy tensor. Sect. HI develops several mathematical results needed to use 
generalized Killing fields in perturbation theory. Some results related to bitensors, 
Riemann normal coordinates, and metric normal tensors used there are collected in the 
appendix. As a simple application of these ideas, the shift in the effective quadrupole 
moment of scalar charge is finally computed in Sect. El 

The metric is chosen to have signature —1,1,1,1. Units are used such that 
G — c — 1, and sign conventions for the curvature follow those of Wald [T3j. Lower 
case Latin letters a, 6, . . . are used to denote abstract indices. Upper case Latin letters 
are used as coordinate indices. A, B, . . . run from to 3, while I, J, . . . can take the 
values 1 to 3. 

2. Mechanics in curved spacetime 

Useful notions of linear and angular momenta are not easily defined for extended bodies 
in generic curved spacetimes. Various concepts have been advanced over the years in 
order to solve different problems. A particularly interesting example was proposed by 
Dixon in the 1970's [TOl jTTllT^l jUl (THj . His momenta arose as part of a general theory 
of multipole moments describing stress-energy tensors T"'' = T^"^' potentially coupled 
to an electromagnetic field Fab- It was shown that - quite surprisingly - all degrees 
of freedom constrained by stress-energy conservation could be encoded in particular 
choices for the body's linear and angular momenta. The ordinary differential equations 
describing the evolution of these two quantities were shown to be equivalent to the 
partial differential equation 

V^T,fc = FabJ\ (2) 

where J° is the body's conserved electromagnetic current density. The time 
evolution of the quadrupole and higher moments is not constrained by stress-energy 
conservatiorQ. Knowing these moments together with the linear and angular momenta 

I There are additional subtleties arising from the assumption that the moments describe a smooth 
stress-energy tensor. The form of the resulting constraints is not completely understood, although 
the analogous problem has been solved for the reduced multipole moments of J" |17| . 
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is sufficient to completely reconstruct a generic stress-energy tensor. 

Both exact integrals and approximate multipole series were derived for the force 
and torque acting on an arbitrary charged object with finite radius. A center of mass 
worldline could also be defined from the momenta. The resulting laws of motion 
have had an important influence on the understanding of test body motion in general 
relativity. The formalism also applies to objects with strong self-fields, although it is 
then much more difficult to approximate the motion in a useful way. 

Despite this, the consequences of Dixon's definitions have been successfully 
studied for a wide class of self-interacting extended charges in flat spacetime [T^ . The 
resulting electromagnetic self-force agreed with expectations only under very strong 
assumptions. It was later understood that the problematic terms were related to 
effective renormalizations of all components of an object's linear and angular momenta 
[6]. This is much more difficult to recognize "by inspection" than the classical mass 
shifts appearing in simple cases like those described by H]). 

Dixon's formalism has recently been modified in order to much more easily and 
naturally take into account these kinds of effects. All momentum renormalizations are 
included at the outset using arguments that do not involve any perturbative expansions 
of the self-field. Scalar and electromagnetic versions of this formalism have both been 
developed in fixed spacetimes [6l [7] . Similar ideas have also been applied in general 
relativity linearized off of an arbitrary vacuum background [8 . Only the results of the 
scalar case will be reviewed here. 

2.1. Generalized symmetries 

If exact Killing fields do not exist, it becomes convenient to define momenta associated 
with vector fields that are "nearly Killing" near a given worldline. Such objects require 
the specification of a frame; namely a timelike worldline T together with a set of 
spacelike hypersurfaces {Ss|Vs G M} foliating a worldtube W D T. Any point in a 
given Es is assumed to be connected to 7.; = Fn Ss by a unique geodesic lying entirely 
within that hypersurface. 

Given such a frame, a set of generalized Killing fields (GKFs) may be defined 
within W. These share several properties with ordinary Killing fields. Most 
importantly, knowledge of a GKF and its (antisymmetric) first derivative at any point 
on F fixes that vector field throughout W. The dependence on these initial data is 
linear and non-degenerate. For a given frame, there therefore exist exactly 10 linearly 
independent GKFs in every four dimensional spacetime. They satisfy 

^idablr = ^aC^gtclr = 0. (3) 

Any ordinary Killing fields that may exist are also GKFs. 

These were the only properties used in [6l [7] in order to derive laws of motion for 
matter distributions coupled to their own scalar and electromagnetic fields. They do 
not provide a unique definition, however. In order to explicitly evaluate the self-force 
beyond its lowest order, the GKFs proposed in [9] will be adopted hercd. They are 
derived from a one-parameter family of Jacobi fields 4'°'{x,s). Associated with each 
point 7s on F are two tensors Aais) and Bab = S[o6](s) related to the Jacobi fields via 

A(s) = i>ahs,s), Babis) = Vatpbils.s). (4) 

§ In [9], a larger class of vector fields were introduced as generalized afliine collineations (GACs). The 
objects needed here are special cases originally referred to as Killing-type GACs. Following [6l[7], we 
simplify the terminology and relabel these as GKFs. 
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Using these quantities as initial data, ip°'{x,s) may be found away from 7^ by 
integrating 



Rabc'^fr^d = (5) 



along an affinely-parameterized geodesic y{u) connecting that point to x. 

Eq. ([5]) is just the geodesic deviation or Jacobi equation. It is easily verified to 
be equivalent to the two first-order equations 

DV-a dy"- dy-(BBab^^ d^y^A . fa^ 

—\ -^ba = 0, -3— — h Rabc -T-Vd = 0, (6) 

du du du \ du du I 

together with initial conditions V'a(7s,s) — Aais) and Babijs,s) = Sah(s). Both of 
these forms have simple interpretations. The first-order equations are very nearly 
the Killing transport equations that relate genuine Killing fields at different points. 
The only difference is the overall factor of y"" in the second half of (O. Without this 
projection, the full Killing transport equations would lead to inconsistently propagated 
derivatives. 

Next, let t{x) be a "time function" that returns which leaf of the foliation its 
argument belongs to. This means that t{x) — s for every a; G Ss- A GKF may 
now be derived from ip"" by setting 

Cix)=r{x,T{x)). (7) 

It may be shown that ^ implies 

^-7Xa = 0, ^ + R^,/^^^A, = 0. (8) 
ds ds 

These are the full Killing transport equations, and should be contrasted with 

They specify a unique family of Jacobi fields - and therefore a unique GKF - for any 

choice of Aaiso) and Bab = B[ab]{so) given at some point € F. 

To summarize, a GKF may be obtained at any point in W from initial data Aa{so) 

and Sah(so). This is done by first integrating ([5]) along an appropriate segment of F, 

and then ^ out to the point of interest. The resulting objects satisfy all of the 

properties discussed above. They are also related to the initial data via 

Uls) = Aa{s), Va6(7.) = Bab{s) = 6[,b](s). (9) 

Other consequences of this definition are discussed in [9] . 



2.2. Momentum 

Given some extended matter distribution with a spatially-compact worldtube W C W, 
we represent its total momentum with respect to the given frame by a linear functional 
V^is) of the GKFs at each point on F. If, for example, Z"" = d/dZ were chosen to be 
a translational Killing field in some spacetime, Vz would represent the Z-component 
of the body's linear momentum. It is more common to represent linear and angular 
momenta by tensor fields Pq(s) and Sab = '5'[a(,](s) on F. These objects are related to 
P{ via 

=Pa(s)r(7.) + Is-'maUls). (10) 

Working with the functional on the left-hand side provides a number of simplifications, 
although it is completely equivalent to using pa and Sab- 
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Useful definitions for these objects have been derived in 6, and [2, for matter 
coupled to scalar and electromagnetic fields, respectively. For simplicity, only massless, 
minimally-coupled scalar fields are considered here. The momenta then have the fornljj 

n - / T\edSa+S^, (11) 

where Tab is the stress-energy of the body in question. £^ represents a contribution 
from the body's S-type self-field. It is defined by (|22p below. 

If the only source of energy or momentum immediately outside of the body is the 
scalar field cf), the stress-energy tensor appearing in pip is given by 

Tab = T'af - ^ (V,0Vb<^ - i5a6V,</)V^0) (12) 

in a neighborhood of W. The total stress-energy tensor T*^* is conserved: 

If p is the body's scalar charge density, its self-field is defined to be the retarded 
solution to 

V'Va^ = -47rp. (14) 
The remaining field is labeled "external;" i.e. 

= 0cxt + 0rct = 0oxt + (/"sclt- (15) 

This implies that 0cxt is a solution to the homogeneous version of lfT4|) in a 
neighborhood of the body's worldtube. The scalar charge density p is assumed to 
have similar support properties to Tab- More precisely, 

W^W = supp Tab 3 supp p. (16) 

The self-field can be further split using a Detweiler- Whiting decomposition [SJ [TS] . 
Define S and R-fields satisfying 

4>i-ct = 4>S + 0R- (17) 

The two quantities on the right-hand side are referred to as the singular (S) and regular 
(R) components of (pseif: respectively. Despite their names, both are mathematically 
well-behaved for physically reasonable charge distributions. Roughly speaking, (jjg 
represents the "Coulomb-like" bound field and </)r the radiation. The latter object 
satisfies the homogeneous wave equation, and is usually responsible for most of the 
non-inertial self-force. These objects are defined via Green functions GR(x,a;') and 
Gs{x,x') using expressions like 

Mx) = / p{x')GKix,x')dV'. (18) 
Jw 

They satisfy 

VVaGR = 0, V"V,Gs - -4n6{x, x'), (19) 

and 

Gsix,x')^Gs{x,x). (20) 
The S-type Green function also vanishes if its arguments are timelike-separated. 

II The notation here has changed somewhat from (6][7]. There, the first term of lllll l by itself was 
referred to as . The sum was denoted by . 
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Two more brief definitions are necessary before writing down the self-field's 
contribution to a body's momentum. Bisect the body's worldtube with E^. Let 
the halves to the future and past of this hypersurface be denoted by E+ and Ej, 
respectively. Also define the field "due to" charge in a spacetime region A by 

0s[A](:e)- / p{x')Gs{x,x')dV'. (2f) 

J A 

The scalar component of the self-momentum appearing in ()f ip is then [6] 

= I dVpC^M^;] - dVpC^MK]^ • (22) 

Despite the unbounded volumes of E+ and E^T , the support of Gs restricts the integrals 
here to finite times into the past and future of order the body's diameter. These 
definitions reduce to standard results in stationary systems [6l[7]- Regardless, the 
total momentum V^^ in (fTTj) has now been fixed in terms of p, gab, and Tab- 

2. 3. Laws of motion 

A body's bulk motion can be interpreted to mean the behavior of its momenta and an 
associated center of mass worldline. We therefore seek a method for determining 7^{(s) 
at arbitrary times given appropriate initial data. The relevant differential equations 
follow from stress-energy conservation and the scalar wave equation p4|) . 

Once they are known, rates of change for the ordinary momenta Pa and Sab are 
easily extracted. Differentiating (fTO|) and applying fSj shows that 

= {pa + lRabcdi'sS'"')C + liS"' - 2p["7.'l)V.6- (23) 

It is clear from this that measures the degree to which the evolution of the body's 
momenta deviates from the Papapetrou equations describing the motion of a spinning 
(but otherwise structureless) test particle. For a fixed GKF it may be thought of 
as a particular linear combination of the net force and torque. 

is derived in [5]. It will be useful here to split it into three parts. Let 

= .Ff + + TI . (24) 

The first term on the right-hand side represents the "gravitational force" (or torque). 
If t° is the time evolution vector field for the foliation {E^}, it is given by 

^r"(s) = i ^ T'^'^{x)C^gbc{xnx)dSa. (25) 

This couples to the quadrupole and higher moments of T"^, and is a generalization of 
effects well-known even in Newtonian gravity. 

Next, J^^°'^ denotes the ordinary force due to the scalar field ^hom = ^cxt + (f'R- 

T^°^{s) = I p{x)C^(j)^,,-^{x)t'^{x)dSa. (26) 

Note that 0hom satisfies the homogeneous wave equation V°Va'/>hom = in a 
neighborhood of W . In most cases, it therefore varies slowly inside the body. It 
is then possible to expand (|26|) in a multipole series about 7^. To lowest order, this 
results in T^°"^{s) ~ (Z>C40hom(7s)- Standard results for the scalar self-force are then 
straightforward to obtain [6]. 
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Lastly, denotes the force exerted by the (inhomogeneous) S-type self-field. It 
is usually very small, as the vast majority of such effects have already been absorbed 
into the definition of . It was shown in ^ that 

= i / cl5X(x) / dV'p{x)p{x')C^Gs{x,x'). (27) 

Lie derivatives of two-point tensors are understood to act independently on each 
argument, so C^Gs = £,°''^aGs -l-C ^a'Gs- The remainder of this paper will focus on 
the consequences of (P7| and its similarity to the gravitational force T^"'^^ . 

This review has only discussed the evolution of a body's linear and angular 
momenta. The choice of frame {F, Sg} has been left open. One useful possibility 
involves introducing center of mass conditions [121 119| . The formalism presented here 
then leads directly to self-contained laws of motion in a more traditional sense. The 
details of such a process will not be needed below. The required steps are described 
in detail in various sources (see, e.g., [71 fT2l [T9] ) . 



3. Effective stress-energy tensors 

It follows from (|27p that the self-force and self-torque due to the S field are determined 
by C^Gg. The singular Detweiler- Whiting Green function is a purely geometric 
construction, so its Lie derivatives with respect to must be linear functionals of 
C^gab- This is also true of Tf. It may be seen more directly by using the field 
equation to replace p{x) in ((27|) by — V''Va0s/47r, and then integrating by parts. 
This requires two intermediate identities. First, note that the second derivatives of 
any vector field can be written in the form [3] 

VfeVa^c = -RcofcdC* + ^(a^i9b)c — -j^ cC.(,9ab- (28) 

This generalizes a result that is well-known for Killing fields (in which case the last 
two terms here vanish). It may be used to obtain a wave equation for C^Gs- Directly 
differentiating lfT9|l shows that 

W^WaC^Gs{x,x') = [W^W'Gs + 27T6{x,x')g-'']£^gab 

+ [V'C^gab - ^Va(ff'^£jg6c)]V'^Gs. (29) 

Combining this with (P7|) gives 

+ / IldSa + f Ilt'dSab, (30) 

where 

^« f drp'[(V>s/:«Gs - ^sV-^^Gs) 

+ ^sCg^'V^Gs - ^g''^W''Gs)C^gbc]. (31) 

The first line of pop has a very simple interpretation. The quantity in parentheses is 
the stress-energy tensor of (j)s (multiplied by An). Comparison with (|^ shows that 
there is a direct sense in which the stress-energy tensor of the body's self-field adds to 
that of its material components. 
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This is especially true in static, asymptotically flat systems. Temporarily suppose 
that there exists a timelike Killing field T", and that each T,g is orthogonal to it. 
Also set CtP = ^t4's = and 7" cx T°. The time derivative in ([50]) clearly 
vanishes in this case. At large radii r, (/>s decays like l/r. Similarly, there exist 
coordinate systems where the metric looks approximately Minkowski near spatial 
infinity: gAB ^ Vab + 0(l/r). Killing fields in flat spacetime can grow like r^. 
If the GKFs can be defined at large distances - which is by no means guaranteed 
- they might grow at similar rates. Assuming this to be true, C^gAB ~ 0(r°) and 
^ 0(f/r^). It follows that the integral over d'Es in ([50)1 will vanish if Ss is pushed 
to spatial infinity. The only remaining term is the volume integral. It implies that the 
body responds to the background gravitational field as though it had a stress-energy 
tensor 

j.ab ^T-o + l- (^V^^sVVs - ^.g^'V^^sV.^s) . (32) 

The extra terms here are exactly those expected for the stress-energy tensor of an 
(isolated) scalar field 0s- 

While suggestive, there are several problems with this result. First, it was only 
derived for static systems. It is not particularly clear that ((32)) can give the body's 
effective stress-energy tensor as it would be inferred by observing responses to different 
spacetime curvatures. As already alluded to, the derivation even in the static case 
is suspect due to both Gs and the GKFs generically being undefined at very large 
distances. Further conceptual and technical problems also arise from the infinite 
extent of the stress-energy tensor associated with (j)s ■ It becomes very difficult to define 
effective multipole moments, for example. This is especially true if the spacetime is 
not Ricci-fiat. Many of these problems can be at least partially circumvented at low 
orders in perturbation theory by using scaling arguments like those in [20] . It appears 
unlikely, however, that this kind of reasoning can be extended very far. No such 
attempt will be made here. 

In general, either (|77|) or (|5D|) (with boundary terms) may be used to rigorously 
compute J-^. The latter equation was useful to motivate results like (15^ . although 
it is simpler to use (j27p in calculations. The crucial point is to note - as already 
mentioned - that C^Gs{x, x') is linearly dependent on C^gab in some region near x 
and x' . If the two arguments of this Green function lie sufficiently close together, it is 
well-known to take on the Hadamard form [18 

Gsix,x')^^[A'/'S{a)-Veia)]. (33) 

6 and O are the usual Dirac and Heaviside distributions, respectively. Synge's world 
function a{x,x') = a{x',x) is defined in the appendix via (jA.l|) . The scalarized van 
Vleck determinant is 

A{x,x') = Aix\x) = -^^t^^ (34) 

where we have used the standard notation aaa' = VaVa'cr. No similarly simple formula 
is known for the tail V{x, x') = V{x', x). It is, however, a smooth biscalar determined 
entirely by the spacetime geometry near x and x' . 

Lie derivatives of a, A and V with respect to are all linear functionals of C^gab- 
It immediately follows from (|A.2p that 



C^a{x,x) = - / du- —C 

2 Jq au du 



e.9ab{y), (35) 
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for example. As in (|A.ip . y{u) is an affinely parameterized geodesic connecting x and 
x' . Also note that a direct differentiation of ([51)1 gives 

\/alnA = -H'''ba''ab', (36) 

where 

H'''a{x,x') = [-a\,{x,x')]-\ (37) 

The "—1" in this last equation denotes a matrix inverse, which is always assumed to 
exist in the regions of interest. It immediately follows from ([36]) that 

In A(x, x') = -[Vae + Va'r' + i/"'"VaVa'£5CT(a;, x')]. (38) 
2Va^° = g'^^C^gab, so, with the help of ([55)) . iZ^A can also be written entirely in terms 
of C^^gab- 

Such explicit demonstrations do not seem to be easily obtained for V{x,x'). 
Regardless, the tail is a smooth solution to the homogeneous wave equation: V^V aV = 
0. It can be uniquely specified via purely geometric boundary data given on a null 
cone with vertex x or x' |18j . It is therefore conceptually clear that V is an entirely 
geometric quantity whose Lie derivatives may also be written in terms of C^gab- This 
is supported by direct Taylor expansions of V that have been computed to fairly high 
order in the literature [HJ [55] . 

It is now clear from ([27p that J-^{s) can be written as a linear functional of C^gab 
as it appears in a certain finite neighborhood of Es n W. Inspection of (|25p shows 
that the gravitational force J^^^^^ (s) is also linear in C^gab, although only on Eg H W. 
Despite this difference, the two forces may still be related if the metric is analytic 
(in a Riemann normal coordinate system with origin 7s). In this case, C^gab can be 
expanded in a Taylor series about jg. As will be discussed in more detail in Sect. [4] 
below, every term will involve a quantity C^gab,ci---c„{ls) for some n>2. gab,ci---c,i is 
the n**^ metric normal tensor (or tensor extension of gab)^ f^nd is defined by (|A.8[1 . 

Using the resulting expansion in ([55)1 or ([57]) yields infinite series for and 
J-"!, respectively. A term involving C^gab,ci---c„{7s) clearly must be contracted with a 
tensor of rank n + 2. In the case of the gravitational force, this tensor is the 2"-pole 
moment /ci - c„abj-g-j of T"'' (divided by 2-n\). See Sect. l5.1l below. Similar coefficients 
in the series for J^^ may be identified as effective shifts in these multipole moments. 
This means that 

1 °° 1 

^r^(s) + :fI{s) -,i:^-'-^\s)C^9ab,c,...,Sls) (39) 

n=2 

for some effective multipole moments I^^""""'' = ri---'^r,ab_^_gjc,-c„ab^ rpj^jg should be 
compared with the purely gravitational force ([55]) . Formulae for the "bare" moments 
JCl■■■c,^ab jjj^ terms of the matter fields may be found in the literature [TT], and are 
rederived here in Sect. 15.11 The shifts (J/ci - Cnafi (j^g ^jjg self-field, and may be 
calculated from ([27]) . Much of the remainder of this paper will be concerned with this 
process. 

It should be noted that similar results are obtained if the metric is not analytic, 
but still varies slowly near Eg n W. In this case, it is often possible to write down 
asymptotic series for the gravitational force and torque that are very similar to those 
just described [14]. The only difference is that the series becomes an approximation 
that must be cut off at finite order. The same comments also apply to 
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If the metric is not smooth near the body, there needn't be any relation between 
^l'^'^^ and Consider, for example, a spacetime that is everywhere flat except on a 
null hypersurface. As the gravitational wave approaches the body in question, there 
will be a very brief period during which the gravitational force exactly vanishes, yet 
^1 ^ 0. Momentum changes or center of mass motions predicted to arise on such short 
timescales are unlikely to have any physical relevance. This paper will therefore be 
concerned only with systems where C^gat can be approximated by high order Taylor 
series in all regions of interest. The requisite smoothness assumptions are very similar 
to those used in the ordinary multipole expansion as it relates to test body motion. 
That problem has been discussed in detail by Dixon [M]. 



4. Expanding C^gab 



A strong case has now been made that the singular scalar self-field mimics the ordinary 
effects of curvature felt even by extended test masses. We now explicitly expand 
^(,gab{x) as outlined above. This is fundamental to calculating the shifts in a body's 
multipole moments. 

Many of the equations used in this section involve multiple spacetime points. 
Indicating which indices are associated with which points is notationally cumbersome. 
A more consistent use of primes will therefore be made in order to avoid confusion. 
Simpler equations will have primes removed when their meanings are clear. Points on 
r will be denoted by symbols that are left unprimed in most cases. 

If one is interested in finding a generalized force 7^^(s) at time s, it is useful to 
expand Lie derivatives of the metric about the point 7^ e F. A series of this sort for 
gab is given by (|A.9p . It involves "separation vectors" = —a°'{x,x') and metric 
normal tensors gab,ci - cn- Lie derivatives may be computed directly. First note that 



C^ga'b'ix') = 



(40) 



^ipga'b' 

{x'.r{x')) 

which follows from ([7]). V'" is the Jacobi field from which (x') is to be 

constructed. It is implicit that the second argument in this field is to be evaluated at 
t(x') only after all differentiations have been performed. This means that 

for example. 

The first term in (|40|) is easily evaluated using (IA.9P and the identity [9] 

C^a^{x',js) = C^<y\'ix',is) = 0. (42) 

This notation is somewhat ambiguous regarding the arguments of the Jacobi fields. 
Fully expanding the left-hand side, 

C^a'^ix', js) = i'"' {x', s)a\, {x' , 7,) + s)a\{x' , 7,) 

- Vf,V"(7., V(a;',7.)- 
This is valid for any s. Combining these equations with ([3]) yields 

^ 1 

f^a'CT fa' / , —X'^^ ■ ■ ■ X'^"£cgab,ci - 



C^ga'b'{x') 



n=2 



(43) 



(44) 



J {x 



From the perspective of this is only needed in the special case s = t{x'). 
valid for general s, however. 



It is 
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Evaluating dtpa' /ds in (I40p requires some additional care. A direct differentiation 
of with respect to s while using ^ and ([5]) shows that 



(45) 



where 7?" „ is defined by ([57]) . Also note that 

a\ = a''a\c' + (r\,a''b, (46) 

which follows from (|A.13|) . If y"{u) is an affinely parameterized geodesic connecting 
7 to a;', both of these results together with (|42| imply the "transport equation" 

d 







(47) 

(a"(«),7s) 

Normalizing u such that y"{Q) — js and ?/"(l) — x' , an integral of this equation yields 



xa/'{y",js)C^9d"f"{y"). (48) 
As with (UH), this applies for any s (as long as 7s is not too far from x'). Making use 
of that equation here yields a series for dV'° /ds in terms of the metric normal tensors 
at 7s. 

Using the result together with (|40| and (|44|) finally gives 



A5a'!,'(a^') - 



(") 



{x',t{x')) 



where 



(n) 

e 



abed/ 



ix',j,) = {n-l) 



duu 



n~2 



—af „{c „d) „hh 



(j/"(m),7) 



for all n > 2. In flat spacetime, 



(n) 

e 



9 " 9 



(49) 



(50) 



(51) 



These expressions succeed in expanding C^9ab{x) up to A^th order about ^t(x)- They 
generalize and unify results found (by direct computation) in ^9' . They are all that is 
needed to expand J-^^'^"^ in a multipole series. 

Unfortunately, computing JF^ requires knowledge of C^gab{x) at points x not 
contained in T^T^^^y Simple closed-form expressions do not seem easy to obtain in such 
cases. Despite this, it is still possible to find a useful approximation. Equation (|49|) is 
of the form 

C^9a'b'{x') ^ E - ^-'b' (^'.^(2^0), (52) 



where 



(") 

Q a'b' ix',s) 



n — I 



(53) 



{x',s} 
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Note that the arguments here are {x',s) rather than the {x',t{x')) appearing in 
Each of these coefficients may be further expanded as 

g ,,,, [x', r{x')) J2 ^ ^ ^, ' Q ab ix', s) (54) 
e=o 

for some L > 0. This wiU require L derivatives of C^gab,ci---c„{'^s)- Using ([5]), 

^^igab,ci - c„{ls) ^ itC^y dgab,c^-cAls)- (55) 

The Thomas replacement theorem discussed in the appendix impHes that ^ dgab,cx - cn 
can be written in terms of (undifferentiated) metric normal tensors. It foUows that 
s-derivatives of C^gab,ci - c„ are always proportional to Lie derivatives of other metric 
normal tensors. Carrying out all of the expansions here therefore leads to a series for 
t^(,9a'b' in which each term is proportional to some C^gab,cx---Cn{ns)- 

There does not appear to be any simple way to compute this exactly. It is, 

(") 

however, possible to obtain an approximation for Q a'b' ip' ,t{x')) that's valid to n 
order in powers of the separation between x' and 7s . Suppose that each hypersurface 
Eg is formed by combining all geodesies passing through 7s and orthogonal to a unit 
timelike vector n°(s). If this satisfies na7s = ~li have 

Vqt(7s) = -na{s). (56) 

Differentiating (j54[) and using standard results for the coincidence limits of various 
bitensors, it is now possible to show that 



(n) 

g^'b' {x',t{x')) 



„(a c) I rb _^_A,b^ I yd-i yd-n 

a'O 'b' \o c - „ _ Y ^s^'' ) ^± ' ■ ■ "^-L 

X ^^9ab,di - d„ 



+ 0(X"+i). (57) 

(x',s) 

This only depends on the "projected separation" 

Xi(x',7s) - -[5\ + j:nb{s)]a\x',js) (58) 

from the origin 7s. In conjunction with ((52)) . it will be fundamental to estimating the 
effective shift in the quadrupole (n = 2) moment of T"''. 

5. Multipole expansions 

5.1. Gravitational force and torque 

As claimed above, it follows from (PS)) and that the gravitational force and torque 
may be expanded in a series of the form 

1 ^ 1 

^f-(s) ^ - ^ -/-■--''(,)/:^g,,^,^...,J^,). (59) 

n=2 

jci---c„ab^g^ represents the 2^-pole moment of the body's stress-energy tensor at time 
s. Given (lA.lOp and (jA.Up . it must satisfy 

jci---c„ab _ j{ci---c„)ab _ jci---c„(ab) ^gQ-j 

and 

j{ci---Cna)b _ jci(c2---c„ab) _ q (6^) 
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Following [m [13] , it is convenient to define an alternative set of moments 

jfl---Uabcd _ jfi---f„[a[cb]d] ^g2) 

The notation on the right-hand side indicates independent (binary) antisymmetriza- 
tions with respect to the index pairs (a, b) and (c, d). It immediately follows from this 
definition that 

jfl---fnabcd _ jfi---frilab]cd _ jfi---fnab[cd] _ j{fi---fn)abcd ^gg-j 

Using the properties of /ci - cnofc^ -^g g^jgQ have 

jfl---fnO-[bcd] _ jfi---f„-i[f„ab]cd _ q j-g^-j 

This allows (|62]) to be inverted. Explicitly, [III [H] 



\n+lj 

for all n > 2. 

Using these relations together with ([^^ and ([55)1 gives an explicit formula relating 
the multipole moments to T'^^: 



jh-Uabcd^^^ = [ ASh't^'T'''''' Xf^ ■ --Xf-X^^X 

JT.S 



(66) 

This result has also been derived in llj using different methods. It is easily verified 
that ([^5]) and ([M)) are satisfied. There is an additional identity 

^^^jh-f„abcd ^ ^^jh-U[alcb]d] ^ Q (g7) 

that also holds for all n > 1. 



5. 2. Self-force due to the singular field 

(") 

Applying the same analysis to is more complicated. It requires Q ab (2;,t(x)), 
which is not known exactly in terms of metric normal tensors at 7^.. An approximation 
in powers of distance (divided by curvature scale) was necessary in order to obtain 
(|57p . We shall therefore adopt a similar small-size assumption for the matter under 
consideration. It is also useful to suppose that all timescales associated with the body's 
evolution are much larger than its diameter. 

The singular self-force arises from the degree to which Gs fails to remain invariant 
under the action of various GKFs. C^Gs may be computed to lowest nontrivial order 
by making use of its Hadamard form ([55]) . Combining ([5^ . ([57]) . and (|A.16p . C^gab 
is seen to reduce to 

C^ga'b'ix') = ^a^^'a^K'iS'c ~ 2^ln,)Xix{C^Radfb{ls) + 0(^3). (68) 

Substituting this into ([SS]) gives 

C^a{x, y) = -^-{XlYlXlYt + 2j:{X + YfXlYflufiY X^] + 

X {X^'X'' + X(V^) + Y''Y'')[nf{Y - Xyf}C^Rabcd + 0(^5), (69) 
where X°- = —a°-{x,^s) and y = — CT°(y,7s). All indices here are associated with 7^. 
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A similar calculation using ([55)1 gives 

C^lnA{x,y) = ^{Y ^ Xr{Y ~ X)''£^Rab + 0{X^). (70) 

C^V is only needed to zeroth order in It may be found using the well-known 

coincidence limit [18] 

\imV{x,y)^^R{x). (71) 

An application of Synge's rule (see, e.g., [18]) to this equation shows that 

y) = ^^?^(7.) + 0{X). (72) 

Together, jSHl), ((Ml), COl), and j72]) specify £^03 to lowest nontrivial order. 

Substituting the result into (|27p and performing the relevant integrations becomes 
significantly simpler if a more precise approximation is adopted for the time evolution 
of the matter distribution. Fix the orthonormal tetrad e\ introduced in the appendix 
such that 65(7^) = n°'{s). Also suppose that the remaining spatial triad is Fermi- 
Walker transported along F. For all 7=1,2, 3, let 

£e? + 2n[°n^l5fc,eJ^ =0. (73) 

This allows the introduction of a scaled "radial separation" 

X'{X,x) = -ei(T(a;))a"(a;,7,(,))/A. (74) 

The parameter A > is introduced for later convenience. A time coordinate for a point 
X may be identified as t{x). The four numbers r and XX ^ then define "pseudo- Fermi" 
normal coordinates with respect to the frame {F, Ssfl 

Now consider a 1-parameter family of matter distributions {p, T"''} described by 
the parameter A. Let these satisfy the scaling relations 

p{X,x) ^ \-^"p{X,T,X'), (75) 
T^^(A, x) = X-^^f^^{X, r, A-^), (76) 

in the Fermi normal coordinate system just described, p and are assumed to be 
smooth in all of their arguments, and to have compact support in . a and (3 are real 
constants. The given relations state that the stress-energy tensor and charge density 
both scale down to F in a self-similar way as A — > 0. It should be stressed that this is an 
assumption on the system's behavior in spacetime. While it is almost always possible 
to choose initial data with some type of spatial self-similarity, time evolution needn't 
preserve it. Most importantly, ([75)) - (|76p imply that the body's internal timescales do 
not scale with their spatial extent. It should be kept in mind that this is demonstrably 
false in many interesting systems. The problem is made much worse for certain values 
of a and (3. A sharp existence theorem would be difficult to obtain, however. 

We shall simply assume that reasonable choices have been made. Note, however, 
that the body's proper diameter D scales like A as A ^ 0. Similarly, its scalar charge 
will be proportional to A'^^^". The (bare) mass scales like A"^"^^. These estimates 

^ An ordinary Fermi normal coordinate system uses a foliation orthogonal to its central worldline. A 
more general choice is allowed here. The extra complication is necessary in order to properly define 
a center of mass frame. 
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may be used to reduce the parameter space in several ways. First assume that the 
potentials remain everywhere finite as A ^ 0. This implies that 

a < 1. (77) 

Also suppose that the center of mass acceleration does not diverge in the small-A limit 
we are considering. This means that the charge-to-mass ratio should remain bounded. 
Hence, 



a<f3. 



(78) 



Both of these inequalities together imply that the body's inertia is not predominately 
due to its self-fields. The case a = (3 = 1 was considered in [3 [20]. Except for 
assuming the given inequalities, no particular choice will be made here. 

The singular self- force may now be computed to lowest nontrivial order. 
Substituting the above matter assumptions and the expansion for C^Gg into (|27|) 
gives 



,F|(.) 



12 



1 



n'^n'^WX-yf g'"^\\X-yf 



e'jeix'y' x'^y 
+ 'i^g'^Wx - yf({x - yY{x - y)'' - 5"'\\x - yf) 

X C^Rabcdils) + 0(A«-4"). 



(79) 



The notation || A'jp denotes the usual Euclidean norm J/jA'^A'''. To quadrupole order, 
the gravitational force ([CT]) reduces to 

-J'^^'^'^C^Rabcd- (80) 



6 



jabcd _ Q(^)!>^2f}-^ jg defined via ^ above. Comparison of ^ and ^ shows that 
an object moves as though it had a stress-energy quadrupole moment 

jabcd ^ jabcd ^ gjabcd^ (g^^ 

where 



g jabcd ^ 



A 



7-4a 



x^y^ + y'y 
1 



g'^^'^Wx-yw 



+ 2ff"" \\x-y\\mx- yy {X - 3^)^ - 5"^ \\x - y\\ 



0(A 



8-4aN 



(82) 



It is easily verified that this satisfies the same index symmetries as the "bare" moment 

jahcd. 



^ jabcd ^j[ah\cd ^jab[cd 

^ja[bcd] _ ^j[abc\d _ ^ 



(83) 
(84) 

The physical content of ((82)l is far from clear as it stands. It is useful to split 
the full quadrupole moment into three pieces = 5Q''°'''\ Stt"''"^ = Stt"'^'"'^ and 

ggabcd ^ ^g[ab]cd ^ ^5'a6[cd]_ [Jg] 

g jabcd ^ .s^lajgfcKc^d] _ n["<57r''l"'' - n^^jTr'']"'' + 55"''"''. (85) 
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This decomposition is unique if each tensor here is assumed to be fully orthogonal 
to Stt"''"', and 55'"'"^'* are then interpreted as shifts in the body's mass, 

momentum, and stress quadrupoles, respectively. In general, 

JttI^^'^I = ^S-I^Hrf ^ j^aM ^ 0. (86) 

Considerable simplifications occur if C^Rati'^s) = 0. This will almost always 
be the case in general relativity if - as we have been assuming - the body's own 
metric perturbations are ignored. The introduction of a cosmological constant does 
not change this situation since C^gabi"fs) = by ([3]). Assuming that the Ricci tensor 
can be ignored, traces of the three quadrupole moments decouple from the laws of 
motion. It is therefore useful to define trace-free versions of these tensors satisfying 

gabSQ^F = gabSn^F = gacSS^'p^" = 0. (87) 

The last equation forces SS^'^p'^ to have all of the same algebraic symmetries as a 
three-dimensional Weyl tensor. It therefore vanishes. A "trace-free" version of the 
full quadrupole moment may be defined as 

5 J-^jf^ = -3n['^<5Q^[?n'^l - n^'^S^f - n^'Sn^^'. (88) 

Despite the name, gac^J^'p'^ 7^ 0. It is, however, much simpler than SJ"''"^'^. Forces 
and torques computed using both moments are the same. Letting Cabcd denote the 
Weyl tensor, 

SJ'^^'^'^C^Rabcd — SJ^p'^C^Cabcd (89) 

when C^Rab — 0. 

Equation ^ can be used to show that Jvrf.'j? = 0(A^~''") and 

J d^Xp^se'}e''jiX'X'' - ls''\\Xf) + 0(A«-4"). (90) 

This is exactly the quadrupole moment that might have been guessed from elementary 
considerations together with an assumption that the effective stress-energy tensor was 
given by ([5^ . The same argument could not have been used to deduce the full 
quadrupole moment (|82p needed in non-vacuum spacetimes. In that case, attempting 
to compute full (finite-trace) moments of the field's stress-energy tensor would lead 
to divergent integrals. The approach taken here cannot be bypassed in general. This 
might be an important point when analyzing motion in alternative theories of gravity. 

We have focused on deriving the correction to an object's quadrupole moment. 
The results of Sect. |4] can be used to derive the lowest-order shifts in the octupole 
and higher multipole moments as well. The main obstacle is a lack of knowledge 
regarding the behavior of the tail V{x,x'). Results are available in the literature up 
to fairly high (but finite) order [21] [22] . These could be used to compute shifts in 
the next few multipole moments above the quadrupole. The situation is considerably 
simpler if Rab — 0. In this case, an expansion for the tail will consist of powers 
of X"" contracted against a polynomial in the Weyl tensor and its derivatives. It is 
impossible to construct any such terms that are linear in the curvature, so 1^ ~ O(C^) 
in this case. Such terms would all be higher order than those due to C^a and in 
the approximations adopted here. In vacuum, the tail can therefore be ignored when 
computing the multipole shifts to lowest nontrivial order. 
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6. Discussion 

We have analyzed the contribution of the singular field to the laws of motion describing 
an extended scalar charge in curved spacetime. In most cases, the total force and 
torque is given by 



The multipole series here is asymptotic in general. The maximum value one may 
choose for N depends only on the degree to which the metric may be accurately 
approximated by a Taylor series inside the body. If it is analytic near E^. H in a 
Riemann normal coordinate system with origin 7^, N may be taken as infinite. The 
right-hand side of (pij) then becomes exact (with the usual caveat that any effects of the 
body on the metric have been ignored) . As noted following , the integral involving 
external and R-type scalar fields can usually be expanded in its own multipole series. 
The monopole term is all that is necessary to recover standard results for the self-force 
on a scalar charge [6]. 

Note that all traces of the S field have been absorbed into the effective (or 
renormalized) multipole moments Methods to compute these have been 

discussed in detail in Sects. [3] and [51 An approximation for the quadrupole moment is 
explicitly given by ([5T|l with ([5^ . and ([5^ . The resulting expression is simplified 
significantly if C^Rab = on the body's worldline. In this case, the force due to the S 
field is approximately 



where SQ^p is a symmetric, trace-free tensor given by (j90p . This is exactly what one 
would expect by applying a Newtonian quadrupole definition to the energy density of 
a static S field. No similarly simple interpretation exists when C^Rat 7^ 0, however. 

As discussed in the introduction, it is common to treat the motion of a small body 
as though all observations of its behavior are "external." Masses are then inferred by 
observing the change in acceleration in response to a known change in an external 
force, for example. The same is also true for the spin angular momentum and all 
higher moments of the body's stress-energy tensor. The laws of motion satisfied by an 
extended test particle contain a number of parameters such as these. They may all be 
fit to the observed variations in the motion of a real object as it experiences various 
perturbing influences. The central result of this paper states that once all such fitting 
is complete, there is only one effect that has not been taken into account for a body 
with large self-fields. That is the force and torque exerted by 0r. In a sense, the 
( "interesting" ) self- force and self-torque is given entirely by the regular component of 
the self-field. This is much smaller than the full self-field. The fact that it satisfies a 
vacuum field equation also means that it usually varies much more slowly than 0soif. 
This situation is very similar to the one usually assumed from the start in effective 
field theory treatments of the self-force problem (see, e.g. ^51 124j ). 

It also serves as a direct derivation of a generalized Detweiler- Whiting axiom. 
Our results have proven a slightly modified version of this conjecture to be accurate 
far beyond any regime in which it was originally expected to hold. Up to effects that 
renormalize test particle parameters, there exist linear and angular momenta definable 
for essentially arbitrary extended charges whose instantaneous evolution depends on 




-ci ■ --Cnab 



^^9ab,ci ■••Cn- 



(91) 



(92) 
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the self- field only through The correct momenta are exactly those introduced in 
[HI and reviewed in Sect. [5] above. Recall that the original proposal was intended to 
treat only small "pointlike" objects to lowest order in some approximation scheme fS]. 
In that context, the S field was assumed only to renormalize the body's inertial mass. 
Despite the vast improvement in generality accomplished here, the precise prescription 
for splitting up the self-field has not changed. 

While this paper has focused on the case of a scalar charge moving in a fixed 
curved spacetime, its basic conclusions are easily generalized. The only essential point 
was that the singular self- force depended on C^Gs, where Gs was a Green function 
defined only using geometric objects. In electromagnetism, the analog of (|27l) is 7J 

^1^1 f ^Sat'' f dV'JtJb'C^Gf, (93) 

where Ja is the charge's current density. Gg"* is the singular Detweiler- Whiting Green 
function, which is again defined purely in terms of the geometry. The situation is 
slightly more subtle when studying gravitational self-forces, although the resulting 
expression for the singular self-force is still very similaiQlS]. The arguments presented 
in this paper here are easily extended to all of these cases. 

The calculations carried out in Section [5] may be of observational interest if 
extended to the gravitational case. In that context, future interferometers might 
be able to precisely measure the effects of a neutron star's quadrupole moment on the 
gravitational waves emitted as it spirals into a supermassive black hole. They may also 
be able to do the same for the star's spin. Such measurements could allow one to learn 
something about the star's equation of state. This would require a precise relation 
between the quadrupole moment extracted from the gravitational wave signal and the 
star's stress-energy tensor. It would be interesting to compare the post-Newtonian 
expectations for this to those derived from the present formalism. The nature of the 
"effacement principle" derived here might also be compared to the ones discussed in 
post-Newtonian contexts [STJ [32l [33l [34] . 



Appendix A. Riemann normal coordinates and metric normal tensors 

A spacetime's world function (or Synge's function) u{x,x') — u(x' ,x) is defined to 
be one-half of the squared geodesic distance between its arguments. Let y{u) be the 
unique affinely-parameterized geodesic connecting two nearby points x and x' . If 
2/(0) — X and y{l) — x' , the world function is given by 

^{x,x)^- Au— —gab{v)- (A.l) 

2 Jo du au 

Multiple derivatives of these objects occur frequently, so it is conventional to express 
them by appending indices: VaVt,' VcCr(x, x') is shortened to <Tcb'a{x,x'), for example. 

+ The methods used here depend on the Hnearity of the field equations. The formalism developed 
in |8] therefore treats gravitational self-fields as perturbations satisfying the linearized Einstein 
equation off of some given vacuum background. More general results might be possible by iteratively 
linearizing, or assuming approximate stationarity or axisymmctry. This is not yet known. There is 
also, at present, no formalism that allows both gravitational and scalar (or electromagnetic) self-fields 
to be large simultaneously. 
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First derivatives of a are particularly interesting. They lie tangent to y{u). 
Specifically, 



Y7 t '^ 

\/a(J(X,X ) = 

du 



u=0 



du 



(A.2) 



M=l 



These objects therefore serve as useful "separation vectors." Let e'^{x) (with A ~ 
0, . . . , 3) be an orthonormal tetrad at a particular point x. In a Riemann normal 
coordinate system with origin x, an arbitrary nearby point x' is associated with the 
coordinate values 

X^{x, x') = -e^ix)a''{x, x'). (A.3) 

Differentiating this with respect to x' shows that 

V,,X^^-e^a\,. (A.4) 

The left-hand side here clearly reduces to the identity matrix when evaluated in 
Riemann normal coordinates with origin x. It follows that —e^cr'^a' is a covariant 
way of writing the identity. aS-A ^^^^ corresponds to the identity in Riemann 
normal coordinates if H"^ a is defined by ([57|) . 

The metric in these coordinates can therefore be written as a two-point tensor 

field 

Gab{x,x') = e^(x)es(x)if'' a{x,x')H'' b{x,x')ga'b'{x'). (A. 5) 

There is no obstacle to rewriting this object with arguments {x,X^) rather than 
{x,x'). Set 

GAB{x,X''{x,x')) = GABix,x'). 
The iVth-order "covariant Taylor series" of Gab is then 



(A.6) 



^ 1 



Gab 

The coefficients here define tensor fields 

fjab.ci ■■■Cn {x^ 



d"G 



AB 



e^e?e^" 



dx^i ■ ■ ■ axC" 



9x^1 . . . axC" 



(A.7) 



(A.; 



that are independent of the tetrad. They are referred to as metric normal tensors 
[TTl [25] . The general process of looking at the tensor fields defined by coefficients in a 
Taylor series in Riemann normal coordinates is known as tensor extension. 5afc,ci - c„ 
is therefore the n*'' tensor extension of gab- 

Inverting (|A.7P with the help of (|A.3p . (fST]) . (|A.5p . and (jA.Sp . the ordinary metric 
gab{x) may be shown to have the expansion 

^ 1 

ga'b' ix') 0^ a\,a\, V -X^' ■ ■ ■ X^- gab,c,-c^ [x). (A.9) 

We have defined X" = -a-''{x,x') by analogy to (|A.3p . The right-hand side of (|A.9p 
makes use of an origin x that does not appear on the left. It represents the usual 
freedom in Taylor-expanding a function about an arbitrary point. 

It is possible to deduce a number of general properties of the metric normal 
tensors. It is clear from their definition that 



9ab,c\---Cn — 9{ab),Ci---Cj^ — 9ab,{c^---Cn)- 



(A.IO) 
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Less obviously, they also satisfy 

ffa(6,ci---c„) — 5'(a&,ci---c„_i)c„ = 0. (A. 11) 

In p dimensions, the nth-order metric extension has 
p{n-l){n+p - 1)! 
2(n+l)!(p-2)! 

algebraicaUy independent components [25l [26] . 

To prove (IA.11|) . recaU that the world function satisfies [18] 

a^cja = 2(T. (A.13) 
Differentiating with respect to x shows that a°'a' = c". Hence, 

H"' ao-" ^ -a"' . (A.14) 
Combining this with (|A.5[) . 

X'^{x,x')[GABix,x') ~ tjab] ^ 0. (A.15) 

rjAB = Gab{x,0) = GAsiXjX) is the ordinary Minkowski metric. Repeatedly 
differentiating (|A.15p and evaluating the result at the origin recovers (|A.11[) . 
Incidentally, this equation also provides an interesting link between Riemann normal 
coordinates and the Fock-Schwinger gauge of electrodynamics [27] [28]. 

First derivatives of the metric vanish in Riemann normal coordinates, so gab,c = 0. 
More generally, Dixon has found all metric normal tensors to linear order in the 
curvature [11]. For all n > 2, 

' n — V 



9a,b,ci- 



1) V(,3...,„i?|,|,,,,),, + 0(i?2). (A.16) 



This equation is exact if n = 2 or 3. Nonlinear terms appear at higher orders, however. 
It is extremely tedious to compute them by hand, although several methods have been 
developed to make the process conceptually straightforward [571 [55] • For n = 4, 

6 16 ^ 

3ah,ciC2C3C4 = g V(ciC2-^|a|c3C4)f) + Yg"^<^('^iC2 -^|6|c3C4)d- (A-1'^) 

Going to much higher orders than this is best done with a computer algebra package 

m- 

In general, go6,ci - c„ can always be written as a polynomial in the metric, the 
Riemann tensor, and its first n — 2 covariant derivatives (if n > 2). As indicated by 
(|A.16p . the linear term always involves exactly n — 2 covariant derivatives of Rabc^- 
Its remainder can include pieces with at least two fewer derivatives than this. These 
results may be inverted. Any number of covariant derivatives of the Riemann tensor 
may be written as a finite polynomial of metric normal tensors (including the "n = 
extension" gab and its inverse). This is easily seen by noting that ^ fi - f^Rabc^ can 
always be written in terms of partial derivatives of the metric in an arbitrary coordinate 
system. Specifically, this makes sense in a Riemann normal coordinate system. But 
each of these partial derivatives may be replaced by a metric normal tensor on account 
of (|A.5|) and (jA.SP . This means that almost any covariant quantity locally constructed 
from the metric may be written entirely in terms of its extensions. This result is 
essentially Thomas' replacement theorem [25j . 
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